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I.  Introduction 

Electron  and  proton  transfer  reactions  are  among  the  most  fundamental  of  chemical 
processes.  The  important  role  of  the  solvent  activation  in  the  electron  transfer  process  was 
recognized  by  Marcus  many  years  ago. In  particular,  the  concept  of  a  solvent  reaction 
coordinate  and  associated  free  energy  functions  have  played  an  important  role  in  many 
theoretical  and  experimental  studies.  These  free  energy  functions  were  first  calculated  by 
Marcus  using  a  continuum  dielectric  approach^’^  and  statistical  mechanics^  for  electron 
transfer  in  homogeneous  solutions.^  for  electron  transfer  at  metal-liquid"*'^  interfaces,  and 
later  for  semiconductor-liquid  and  liquid-liquid  interfaces.®’’^  This  perspective  on  the  sol¬ 
vent  participation  in  electron  transfer  was  further  developed  by  Levich®  and  Dogonadze.*^ 

Free  energy  functions  for  electron  transfer  reactions  have  been  studied  extensively  us¬ 
ing  Monte  Carlo  and  molecular  dynamics  simulation  techniques.  For  example.  Warshel 
and  co-workers*®  have  studied  the  behavior  of  the  Marcus  free  energy  curves,  particu¬ 
larly  in  the  “inverted"  region.  Kuharski  et  a/.**  performed  extensive  simulations  on  the 
Fe'^^ /Fe^^  redox  pair  in  water  and  provided  conclusive  evidence  that  Marcus  behavior 
would  be  valid  in  realistic  condensed  phase  systems.  Carter  and  Hynes*^  also  performed 
a  molecular  dynamics  study  and  found  larger  cuT'vature  of  the  solvent  free  energy  near  an 
ion  pair  than  near  a  neutral  pair  of  atoms— an  effect  also  investigated  earlier  by  K2ikitani 
and  Mataga.*®  though  both  studies  were  called  into  question  by  Tachiya.’®  More  recently. 
Benjamin*"*  has  calculated  the  free  energy  curves  for  electron  transfer  across  the  inter¬ 
face  between  two  immiscible  liquids,  while  Straus  and  Voth  have  done  the  seune  for  the 
metal-water  interface.*® 

In  the  field  of  proton  transfer,  the  current  understanding  of  the  sol-vent  participation  in 
the  reactive  event  again  has  its  origin  in  the  phenomenology  for  charge  transfer  developed 
by  Marcus.*"’  and  others.®'®  Computational  efforts  have  focused  on  the  calculation  of  the 
solvent  contribution  to  the  activation  free  energy,  studies  of  intramolecular  effects 
eurising  from  the  modulation  of  the  proton  tunneling  by  molecular  vibrations. and 
simulations  of  adiabatic  proton  transfer  in  strongly  hydrogen-bonded  complexes. 

As  opposed  to  simulations,  integral  equations  have  been  employed  to  calculate  sol¬ 
vent  activation  free  energies  for  charge  transfer  reactions  in  only  a  few  cases.  Rossky 


and  co-workers*^  carried  out  extended  Reference  Interaction  Site  Model  (RISM)  integral 
equation*^ calculations  to  examine  the  aqueous  phase  S_\2  reaction  of  a  chloride  ion 
with  methyl  chloride.  Tliis  integral  equation  approach  was  particularly  effective  for  effi¬ 
ciently  examining  the  shape  of  the  free  energy  surface,  the  variation  of  the  results  with 
changes  in  the  solvent  model,  and  the  relative  contribution  of  entropy  to  the  net  free  energy. 
In  the  held  of  electron  transfer.  Kuharski  and  Chandler^^  have  also  applied  the  extended 
RISM  method  to  study  the  solvation  structure  and  thermodynamics  of  di-  and  tri-valent 
iron  ions  in  aqueous  solution.  In  other  integral  equation  studies  of  electron  transfer,  analyt¬ 
ical  expressions  for  the  Born  solvation  free  energy  of  a  pair  of  ions  with  arbitrary  charges, 
unequal  size,  and  hxed  separation  have  been  obtained  via  an  interaction-site  approach  with 
the  mean-spherical  approximation  ^MSA).'*'^  In  the  latter  study,  the  results  for  ions  in  a 
dumbell  solvent,  in  a  dipolar  hard-sphere  solvent,  and  in  an  ionic  solution  in  a  dielectric 
continuum  solvent  were  compared.  The  results  from  the  theory  for  the  elec¬ 

tron  transfer  were  qualitatively  compared  with  simulation  results.  Most  recently,  Fonseca 
pt  have  calculated  the  solvent  force  constant  which  governs  solvation  free  energy 

fluctuations  via  the  reference  hypernetted  chain  (RHNC)  integral  equation  approach^^^ 
for  ionic  solutes  in  model  dipolar-quadrupolar  solvents. 

It  should  be  noted,  however,  that  there  are  numerous  challenges  in  using  existing 
integral  equation  theories  to  calculate  activation  free  energies  for  charge  transfer  reactions. 
For  example,  one  difficulty  is  in  the  appropriate  form  of  the  closure  relation.  In  integral 
equation  theories,  the  closure  enumerates  the  diagrams  which  depict  the  interactions.  Out 
of  all  the  relevant  closures,  the  Percus-Yevick  (PY)  closure  has  proven  most  appropriate 
for  hard-sphere  systems.  Thus,  short  range  interactions  are  reasonably  accounted  for  with 
the  PY  closure.  By  separating  the  bare  potential  into  a  long  range  and  a  short  range 
contribution,  the  MSA  closure  corrects  the  PY  closure  for  weak  attractive  or  repulsive 
tails  to  a  repulsive  core,  but  this  approach  becomes  less  accurate  for  strongly  interacting 
long  range  systems.  For  the  latter  kind  of  system,  the  HNC  closure  has  proven  to  be 
quite  useful  in  RISM  theory.  However,  while  the  HNC  closure  does  introduce  dielectric 
screening  into  the  RISM  theory,  the  dielectric  constant  it  predicts  depends  only  on  the 
intramolecular  correlations  defined  by  the  correlation  function  u.  The  latter  quantity  is 


independent  of  the  intermolecular  correlations  and  thereby  RISM  calculations  may  yield 
an  inaccurate  value  for  the  dielectric  constant. 

Even  though  the  application  of  integral  equation  theories  to  charge  transfer  reactions 
in  polar  solvents  may  be  problematical,  there  is  a  strong  intellectual  and  practical  mo¬ 
tivation  to  do  so.  With  regards  to  the  practical  motivation,  one  needs  only  to  consider 
the  difficulties  inherent  in  computer  simulations  of  systems  with  strong  long  range  forces. 
For  example.  Bader  and  Chandler^^  have  used  computer  simulation  methods  to  study  the 
potential  of  mean  force  between  two  iron  ions  in  a  periodically  replicated  simple  point 
charge  model  (SPC)  of  liquid  water.  Their  numerical  results  were  found  to  depend  sig¬ 
nificantly  on  the  treatment  of  the  long  range  forces.  Specifically,  if  the  long  range  forces 
were  truncated  at  half  the  simulation  box  length,  the  potential  of  mean  force  between  the 
iron  ions  was  found  to  be  nonmonotonic,  exhibiting  an  attraction  at  interionic  distances 
beyond  6  A.  However,  it  was  also  found  by  Bader  and  Chandler  that  if  Ewald  sums  were 
performed,  the  potential  of  mean  force  was  monotonic  and  close  to  a  screened  Coulombic 
type  of  interaction. 

In  the  present  paper,  the  accuracy  of  the  RISM  method  for  calculating  the  solvent 
contribution  to  activation  free  energies  in  electron  and  proton  transfer  processes  is  ex¬ 
plored.  In  the  integral  equation  studies  described  in  the  following  sections,  the  solvent 
system  consists  of  the  water  molecules  and  comprises  the  most  basic  input  for  the  the¬ 
ory.  namely  the  solvent-solvent  correlation  fimctions.  The  solute-solute,  solute-solvent  and 
solvent-solvent  interactions  are  used  as  accurate  inputs  (either  obtained  firom  simulations 
or  analytically  derived)  for  the  appropriate  (in  this  case  HNC)  closure.  The  solvent  con¬ 
tribution  to  the  excess  free  energy  is  then  calculated.  Since  the  ions  are  a  perturbation 
on  the  pure  solvent,  the  quantity  being  calculated,  therefore,  is  the  solvent's  response  to 
the  ions  and  its  corresponding  reorganization.  The  difference  in  excess  free  energy  of  the 
entire  solvent-solute  system  between  its  reactant  state  and  the  transition  state  is  defined 
as  the  solvent  activation  free  energy. 

The  following  sections  of  this  paper  are  organized  as  follows:  In  Section  II,  a  brief 
review  of  the  RISM  equation  is  given.  An  application  of  the  RISM  theory  to  the  model  for 
the  aqueous  Fe^''‘/Fe^'''  electron  transfer  reaction  is  then  given  in  Sec.  III.  Next.  Sec.  IV 


contains  an  application  of  the  RISM  theory  to  a  model  aqueous  proton  transfer  reaction 
which  differs  signihcantly  from  the  electron  transfer  system.  Concluding  remarks  are  given 
in  Sec.  \’. 

11.  RISM  Theory 

.4.  Formalism 

The  RISM  integral  equation.*®'"^^  in  its  extended  form.^^-^®  can  be  applied  to  liquid 
solutions  with  long  range  forces.  For  example.  Pettitt  and  Rossky^^  have  formulated  and 
applied  the  extended  RISM  equation  to  several  three-site  models  of  water.  In  that  study, 
rhc  site-site  cor’''»latir>ri  functions  were  calculated  and  compared  to  available  computer  sim¬ 
ulation  results.  Furthermore,  the  variation  of  liquid  state  structure  with  the  model-site 
charu'‘  was  examined,  and  it  was  found  that  the  integral  equation  approach  provides  a 
correct  qualitative  description  of  the  liquid  structure.  Rossky  and  coworkers^®”^®  have 
also  provided  considerable  evidence  that  the  extended  RISM  approach  provides  a  reason¬ 
able  alternative  to  computer  simulations  for  the  analysis  of  aqueous  solutions  and  ionic 
solvation. 

The  explicit  RISM  equation^®"®®  for  a  molecular  fluid  of  molecular  number  density  p 
is  given  by 

phayP  =  U  *  Cay  *  UJ  +  UJ  *  Cay  *  phayp  (1) 

where  all  matrices  are  site-site  labeled;  a.7  being  two  distinct  sites  within  the  same  species, 
or  different  species,  as  the  case  may  be.  h  is  the  matrix  of  intermolecular  pair  correlation 
functions,  c  is  the  matrix  of  direct  correlation  functions,  p  is  a  diagonal  matrix  of  site 
densities,  and  u  is  an  intramolecular  correlation  function.  The  latter  quantity  is  defined 
as 

OJaMyM'  =  PM^M  M'i^cMyM' 6{raM  -  CyM')  4-  (1  -  SaMyM')  SaMyM' (2) 

for  site  a  in  molecule  M  and  site  7  in  molecule  M'.  Here,  6aMyM'  is  the  Kronecker  delta, 
S(t)  is  a  delta  function,  and  sir)  is  the  intramolecular  pair  distribution  function  between 


distinct  sites  in  the  same  molecule.  For  rigid  bonds,  the  latter  function  is  taken  to  be  a 
delta  function  in  r  space.  In  the  above  notation,  an  analog  of  the  HNC  equation  constitutes 
the  closure  relation  and  is  written  as 

Cq-^y  4“  0cfy  ^or'V  ln[l  “I"  /lof-y]  (3) 

B.  Numerical  Solution  of  the  RISM  equation 

The  set  of  RISM  equations  to  be  solved  now  will  be  discussed.  As  previous  work 
has  shown.^®  the  correlation  matrices  can  be  considered  site-ordered  according  to  different 
classes  of  species.  The  species  in  the  present  case  are  the  solvent  or  the  solute.  Equation 
(1)  then  separates  into  three  equations  for  h"’'.  h"’'.  h““.  which  are  the  solvent-solvent, 
solvent-solute  and  solute-solute  correlation  function  matrices,  respectively.  The  following 
set  of  coupled  integral  equations  are  then  obtained: 

phZp  ^  0J*cZ*uj  +  ui*cZ*phZp  .  (4) 

K'^P  =  +  ■  (5) 

.  (6) 

where  p  refers  only  to  the  solvent  (i.e.,  in  the  present  applications,  the  density  of  the  solute 
is  entirely  negligible  compared  to  the  solvent).  In  Eqs.  (4),  (5)  and  (6),  ha-,  and  Ca->  may 
be  the  solvent-solvent,  solute-solvent,  or  solute-solute  pair  correlations.  Since  each  site 
may  carry  a  charge,  it  is  most  convenient  to  cast  these  equations  in  renormalized  form 
and  then  use  a  closure  which  takes  advantage  of  this  renormalization.  The  details  of  this 
method  was  done  by  Rossky  and  co-workers  in  several  of  their  papers.^^-^®'^®"'’*’  With  the 
closure  relation  in  Eq.  (3),  one  can  solve  for  Eqs.  (4),  (5)  and  (6)  sequentially.  That  is. 
first  Eq.(4)  is  solved  for  the  pure  solvent-solvent  correlations.  After  that,  the  solution  of 
Eq.(4)  is  used  as  input  for  Eq.(5)  and  one  solves  for  the  total  solvent -solute  pair  correlation 
functions.  Next,  the  solution  of  Eq.(5)  is  used  to  get  the  total  solute-solute  pair  correlation 
functions  in  Eq.  (6). 
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The  solvation  free  energy  Usoiv  is  obtained  from  the  Singer-Chandler  relationship.  “ 
It  is  the  solvation  contribution  to  the  solute  chemical  potentials  and  is  given  by 

1 

dr[cas{r)  -  -h'i^ir)] 

1 

+  E  /_  as  [k)x,3'(k)Cs'-yik)eip[ik.T)  (', ) 

where  a  and  7  are  indices  for  the  solute  while  s  and  s'  are  indices  for  the  sites  of  the 
solvent.  The  correlation  function  \s3'  in  Eq.  (7)  is  expressed  within  the  context  of  the 
right -hand-side  of  Eq.  (4)  as  uj  -t-ph^XP-  An  alternative  solvation  free  energy-  expression.  *^ 
which  can  be  derived  from  Gaussian  field  theory,  is  given  by  Eq.  (7)  but  without  the 
( l/2),>i^3(r)  term.  The  latter  expression  may  give  better  results  in  some  cases. but  it  was 
not  found  to  give  substantially  different  results  from  Eq.  (7)  in  the  present  work.  From 
the  last  set  of  pair  correlation  functions  in  Eq.  (6),  one  can  also  determine  the  potential 
of  mean  force  between  the  tw’o  solute  sites  a  and  7.  given  by 


+  .  (8) 

Equations  (4)-(6)  were  aU  expressed  as  matrix  convolution  products  in  real  space.  The 
equations  are  iterative  and  were  expanded  through  fifth  order  in  the  integrals  themselves. 
While  performing  the  calculations,  the  iterative  equations  were  also  expanded  through 
seventh  order.  However,  there  was  not  much  difference  between  the  seventh  order  integrals 
and  fifth  order  expansions.  Each  of  the  three  integral  equations  were  solved  independently 
using  the  appropriate  HNC  closure  self  consistently.  The  system  of  non-linear  equations 
were  numerically  solved  in  real  space  using  standard  iterative  Newton- Raphson  techniques. 
The  numerical  accuracy  for  the  multidimensional  integrals  was  of  the  order  of  10“^  and 
they  were  calculated  by  using  a  256  point  Gaussian  quadrature  subroutine.  The  accuracy 
for  solving  the  roots  of  the  nonfinear  equations  was  of  the  order  of  10~®.  The  number  of 
iterations  needed  for  the  integral  equations  to  converge  ranged  from  45-123.  In  each  case, 
global  convergence  was  achieved  to  a  tolerance  level  of  10“®  for  a  change  of  the  initial 
guess  by  2  orders  of  magnitude. 
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As  an  accuracy  check,  the  same  RISM  code  was  also  used  to  studv  two  systems  unre¬ 
lated  to  the  present  systems.  The  hrst  test  involved  calculating  the  equilibrium  solvation 
tree  energ\"  as  a  function  of  reaction  progress  of  the  aqueous  phase  S\2  ri'action  of  chlo¬ 
ride  with  methyl  chloride.  As  mentioned  earlier  in  the  introduction,  the  same  calculations 
were  performed  earlier  by  Huston.  Rossky  and  Zichi.*^  As  such,  the  solute  and  solvent 
parameters  were  all  obtained  from  the  work  of  Huston  et  al.*^  The  results  of  Ref.  25  were 
reproduced  by  our  RISM  code.  A  second  test  of  the  code  was  to  calculate  structural  prop¬ 
erties  of  XaCl  in  water.  This  work  on  alkali  halides  in  water  was  originallv  performed  by 
Pettitt  and  Rossky.^®  Again.  Eqs.(4)-(6)  were  used  to  hnd  the  pair  correlation  functions 
for  Xa-()  and  Cl-0.  as  well  as  Xa-H  and  Cl-H.  The  model  used  was  the  same  as  in  Ref. 
36.  The  radial  distribution  functions  obtained  from  the  present  code  were  in  agreement 
with  those  of  Ref.  36. 


III.  Application  of  RISM  to  Electron  TVansfer 

.4.  Electron  Transfer  Model 

The  electron  transfer  model  studied  in  this  section  consists  of  the  Fe^^  redox 

pair  in  water.  For  the  most  part,  the  model  is  the  same  as  that  used  by  Kuharski.  ct  a/.‘* 
At  the  reactant  state,  the  charges  of  the  ions  were  fixed  at  -t-2  and  -f3.  At  the  transition 
state  the  charges  on  both  of  the  ions  were  taken  to  be  +2.5.  consistent  with  the  'tight 
binding"  assumption  that  the  electron  density  will  resonate  symmetrically  between  the 
Fe^^  cores.  The  difference  in  the  excess  free  energies  between  the  two  solvated  solute 
charge  distributions  was  taken  to  be  the  solvent  activation  free  energy  in  this  model. 

The  solvent  in  this  model  was  the  simple  point  charge  (SPC)  water  model  of  Berendsen 
et  al.'^‘^  at  a  temperature  of  T  =  298K.  and  an  experimental  number  density  at  1  atm  of 
p  =  0.03334  A“^.  The  total  solvent  pair  potential  in  the  SPC  model  consists  of  a  long 
range  site-site  Coulomb  term  and  a  short  range  site-site  Lennard-Jones  term.  i.e.. 

-/3u- (r)  =  -,5u;7(r)  +  .  (9) 

where  the  pair  potentials  for  water  have  been  decomposed  in  this  notation.  The  short 
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rango  interactions  are  written  as 


The  coehicients  tor 
U'rms  are  -vritten  as 


1  -> 


A2 


and  C^’’  are  listed  in  Table 


I.  In 


(10) 

Eq.  (0).  the  long  range  Coulomb 


o 


L'l’ 

a-y 


(11) 


where  r\  denotes  the  positions  ot  atom  i  of  molecule  a  u  =  O.H.H').  The  partial  charges 
'll!  and  i)ir  are  equal  to  0.41f  and  qa  =  -0.82e  for  ox\'gen.  where  -»  is  the  charge  of  an 


1‘iectron.  To  avoid  an  unphysical  overlap  from  occurring  between  the  hydrogen  and  the 
oxvgen  atoms  in  RIS.  '  calculations,  an  alteration  to  the  SPC  model  has  been  suggested 
l)v  Pcttitt  and  Rossky.^^  This  alteration  amounts  to  adding  a  repulsive  term  eia/r)^'^ 
TO  the  0-H  intermolecular  site-site  potential,  where  f  -  0.8  kcal/mol  and  a-  ~  1,6  A. 
This  modification  is  required  to  ensure  ability  of  the  model  since  without  n  the  total 
water-water  pair  potentials  possess  configurations  with  -x  potential  energv'. 

When  .solvated  in  water,  both  the  ferric  and  ferrous  ions  are  six- fold  coordinated  by 
water  ligands.  The  symmetric  complex  is  octahedral,  with  the  oxygen-ion 

iiond  lengths  of  approximately  2.0  A.  The  oxygen-ion  bond  lengths  arc  approximately  0.13 
A  longer  for  the  symmetric  ferrous  complex  Fe'^^(H20)f5.  An  effective  ferric- water  pair 
interaction  was  constructed  by  Kuharski  et  a/."  that  leads  to  the  same  structural  behavior. 
This  effective  pair  potential  for  the  solute-solvent  interaction  is.  as  usual,  divided  into  a 
short  range  interaction  term 


(12) 


where  Vr,  is  the  position  of  the  iron  ion  and  A  =  6392.7  (kcal/mol) A®,  and  the  long  range 


Coulombic  interaction  of  the  form 


Since  the  valence  ciectron  distribution  in  the  Fe-'  ion  is  locahzed  to  a  great  degree,  as 
noted  earlier,  one  rnav  estimate  the  Fe*  '  fHjOic  interaction  trom  Eq,  \  13)  bv  reducing  the 
charge  trom  3c  to  2c.  That  is.  the  Fe'*  -H-jO  interaction  is  given  by  Eq.  (  13)  with  "3'' 
replaced  by  "2" . 

For  the  solute-solute  {u-u\  interactions,  a  Coulombic  potential  was  used  for  the  in- 
teriomc  terms.  It  is  to  be  noted,  however,  that  the  corresponding  intramolecular  pair 
distribution  function  is  alwavs  a  delta  function  when  the  ions  are  held  rigid.  It  is  also 
important  to  note  that  one  is  rarelv  concerned  with  the  .solute-solute  correlations  per  se. 
but  rather  in  the  solvent -mediated  potential  of  mean  force. 

D  Results  for  Electron  Transfer 

The  RIS.M  result  for  the  activ.nrion  energy'  in  the  Fe^*'  'Fe'^*^  idectron  transfer  problem 
at  a  .separation  of  6.075A  was  found  to  be  2U.8  kcai/mol.  This  result  is  within  591  of  the 
f’xjstmg  simulation  result. '  ‘  The  number  of  oxygen  atoms  in  the  first  solvation  shell  of  each 
iron  ion  was  calculated  to  be  6.2  at  the  t  -ansitlon  state,  which  is  consistent  with  simulation 
result^'  and  known  ligand  coordination.  In  the  reactant  state,  the  number  of  first-shell 
oxygen  atoms  for  the  Fe^*"  and  Fe^*"  ions  was  6.8  and  6.7.  respectively.  Figure  1  .'■hows 
the  Fe^*"  -O  and  Fi'**"  -()  radial  di..trlbutlon  functions  from  the  RISM  calculation,  while 
Fig.  2  shows  the  corresponding  potential  of  mean  force  between  the  two  ions.  Compared 
with  Bader  and  Chandler's  calculations.^'*  Fig.  2  matches  well  with  their  results  obtained 
using  Ewaid  summations.  The  present  results,  however,  do  not  agree  with  earlier  RISM 
calculations  of  the  Fe^'*'/Fe^*‘  solvent  activation  free  energy'.^*  The  reasons  for  this  are 
unclear. 

IV.  Application  of  RISM  to  Proton  Transfer 

A.  Model  Proton  Transfer  System 

This  section  describes  the  application  of  the  extended  RISM  theory  to  a  model  in¬ 
tramolecular  proton  transfer  reaction.  In  principle,  this  problem  is  similar  to  the  electron 
transfer  problem  but  for  two  important  details.  First,  the  ion  distance  is  held  rigid  in  the 
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electron  transfer  case,  whereas  the  proton  transfer  complex  is  llexible.  Thus,  although  the 
basic  tbrmaiism  will  not  change,  the  intramolecular  direct  and  pair  correlation  functions 
of  the  proton  transfer  system  will  be  different.  Second,  the  solute  has  three  sites  with 
charges,  rather  than  onlv  two  as  in  the  electron  transfer  case. 

The  model  proton  transfer  system  consists  of  a  three  atom  solute,  [A  H— A']"*",  where 
A  and  A'  atoms  are  the  same  species,  and  the  complex  is  solvated  in  water. The  three 
RISM  sites  of  the  solute  are  the  atoms  A.  A'  and  the  proton.  The  proton  is  free  to  move 
between  the  A-H-  ■  A'  to  the  A-  •  H— A'  configurations.  The  asymmetric  stretch  mode  is 
taken  to  be  the  reaction  coordinate.  The  solute-solvent,  .solute-solute  and  solvent-solvent 
interactions  are  summarized  below.  The  RISM  equations  in  Eqs  (4)  — (6)  were  then  solved 
with  the  HNC  closure  in  Eq.  (3)  and  the  decomposition  of  the  pair  potentials  in  Eq.  (9). 
The  solvent  activation  free  energy  is  defined  as  the  difference  in  solvation  free  energy  when 
the  proton  is  in  the  reactant  state  to  when  the  proton  is  at  the  transition  state.  As  in 
the  electron  transfer  problem,  the  solvent  contribution  to  the  excess  free  energy  was  found 
from  Eq.  (7). 

The  solvent  in  the  proton  transfer  model  was  the  simple  point  charge  (SPC)  water 
model  at  the  same  temperature  and  number  density  as  before.  For  the  solute-solvent 
interactions,  the  short  range  interactions  between  the  atoms  at  sites  A  (and  A')  and  the 
o.xvgen  atoms  of  water  were  of  the  Lennard- Jones  type,  with  the  coefficients  the  same  as 
fhat  of  ()-()  interactions  in  SPC  water.*"*  Analogous  to  the  extra  O-H  repulsive  term 
added  in  the  solvent -solvent  interactions,  an  additional  repulsive  term  for  the  potential 
energy  between  A  (and  A')  and  the  hydrogen  atoms  of  the  solvent  was  also  added.  The 
proton.  A  and  A'  atoms  interact  with  the  atomic  sites  on  the  water  molecule  by  atom  by 
Coulombic  forces.  The  charge  on  the  solute  sites  was  determined  by  a  charge  swfrching 
procedure  as  a  function  of  the  solute  asymmetric  stretch  coordinate.^**  The  proton  charge 
was  held  at  a  constant  charge  of  -1-0. 5e.  The  reactant  and  product  charges  varied  under 
the  constraint  that  the  entire  solute  system  should  have  a  conserved  charge  of  unity.  The 
charges  of  the  reactant  and  product  sites  were  given  by 

CAiq)  =  ^'>■‘25(1  -  tanhiq/qaw))  ■  (14) 


» 
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eA'(q)  -  0.25(1  +  tanh{q/q^w}) 


(15) 


In  the  above  equations,  e.4  and  e.4'  are  the  charges  of  A  and  A',  respectively,  q  is  the 
reactant  coordinate,  and  q^w  equals  0.125  A. 

UnUke  the  electron  transfer  problem,  the  solute  system  was  not  a  rigid  system.*®  The 
potentials  between  the  protcjiS  aiid  the  A  atoms  were  Morse  oscillators: 

Va-h  =  £)[eip(-2a(lr4Hi  -  b))  -  2[exp{-a{\rAH\  -  5))]  (16) 

In  the  above  equation,  the  parameters  were  chosen  to  be  D  =  21.4  kcal/mol,  a  =  5.4  A~*. 
and  6  =  0.95  A.  The  potential  between  the  A  and  A'  atoms  was  taken  to  be 

Vaa'  =  (l/2)m.u.-L,(r44'  -Te,)'-^  .  (17) 

where  rrir  =  20.0  amu,  ujaa'  =  550  cm~^  and  r^q  —  2.5  A.  A  bending  potential  was  also 
included  to  help  mmntain  the  linear  configuration  of  the  complex.  This  potential  was  of 
the  form 


Vaha'  =  (l/2)fXrul„^  r'Ue  -  Oeqf  •  (18) 

where  =  180®.  pr  —  20.0  amu,  ijJbend  =  1000  cm~*.  and  tq  —  2.5  A. 

B.  Results  for  Proton  Transfer 

The  RISM  result  for  the  solvent  activation  free  energy  in  the  proton  transfer  problem 
was  found  to  be  5.8  kcal/mol,  while  the  simulation  result  is  4.4  kcal/mol.  The  latter  result 
differs  &:om  that  reported  in  Ref.  20  because  432  solvent  molecules  were  used  in  the  present 
case  along  with  an  improved  Monte  Carlo  sampling  algorithm.  In  addition,  the  repulsive 
interaction  term  that  was  added  between  the  interaction  sites  A  and  A'  and  the  hydrogen 
atoms  in  the  RISM  calculation  was  not  included  in  simulations.  Otherwise,  the  simulation 
was  the  same  as  in  Ref.  20. 

The  disagreement  between  the  RISM  and  simulation  results  is  larger  than  in  the 
electron  transfer  case,  though  the  result  is  still  reasonable.  It  should  be  noted  that  the 


rotal  solvent  activation  tree  energj'  is  much  smaller  tor  the  proton  transfer  than  for  the 
Fe*'^  Fe'^'^  electron  transfer^'  or  for  the  Sn'2  reaction'^  (i.e.,  ~  5  versus  ~  20  kcal/mol),  so 
it  IS  much  easier  to  make  a  large  relative  error.  Nevertheless,  the  subtleties  of  describing 
three-body  charge  transfer  processes  with  RISM  theory  remain  to  be  explored. 

V.  Concluding  Remarks 

Extended  RISM  calculations  of  the  solvent  activation  free  energy  for  the  Fe^'''/Fe^''' 
electron  transfer  system  and  a  model  proton  transfer  reaction  have  been  presented.  The 
agreement  with  simulation  results  in  the  electron  transfer  case  was  quite  superior  to  the 
proton  transfer  case.  The  main  differences  between  the  electron  and  proton  transfer  model 
were  the  flexible,  three  site  solute  in  the  latter  case  along  with  a  much  lower  solvent 
activation  free  energy.  In  general,  it  is  perhaps  fair  to  say  that  the  present  RISM  results 
are  mixed,  but  still  rather  encouraging.  Future  research  will  hopefully  shed  more  light  on 
the  use  of  integral  equations  to  predict  solvent  effects  in  condensed  phase  charge  transfer 
reactions. 
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TABLE  I.  Short-Rang  '  Potential  Parameters'*  for  SPC  water. 

—  -  ...  0 

^y0~0)  ^vO-O'l  .-,yO-H)  ^yO-H)  ^yH-H)  ^yH-H) 

Cg  0^2  '-g  Li2  C  12 

-625.731  629624  0.000  225.180  0.000  0.000 

_  • 

“^The  units  for  Cg  are  in  kcal  A®/mol  and  for  C12  are  in 
kcal  A^^/mol. 

i 


Figure  Captions 


Figure  1:  Fernc-oxygen  (solid  line)  and  Ferrous-oxygen  (dashed  line)  radial  distribution  func 
tions  from  the  RISM  calculations. 

Figure  2:  Potential  of  mean  force  for  the  Fe^^/Fe^*"  system  calculated  from  RISM  theory. 
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